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This talk^j is a status report on our study of quantum transport equations relevant 
for baryogenesis computations. Our main finding is that, as a consequence of 
localization in space, the quasiparticle picture of the plasma dynamics breaks down 
at first non-trivial order in gradient expansion. While in this talk we focus on 
bosons, we expect that a similar picture holds for fermions. We then argue that 
the quasiparticle picture is recovered in the adiabatic limit of frequent scattering. 
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1 Introduction and motivation 

A reliable computation of baryon production at the electroweak phase transi- 
tion requires a new formulation of the dynamics of the CP violating sources in 
an out-of-equilibrium relativistic plasma. To this purpose we are working on 
controlled derivation of the relativistic quantum transport equations for weakly 
coupled plasmaLrclll These equations are essential to a systematic treatment 
of the matter-antimatter creation at the electroweak phase transition, since 
they describe the relevant CP violating sources, their transport and dissipa- 
tion in one formalism. Our aim is to apply this formalism to computation of 
the baryon production at the electroweak transition in the Standard Model 
and its supersymmetric extensions, and thus make contact with the upcoming 
particle physics experiments, most notably LHC, DAPHNE, PEP-II, etc. 

An essential ingredient in our method is an expansion in gradients of a 
slowly varying background. It is necessary to go beyond the leading order in 
gradients, since the CP violating sources required by any baryogenesis mech- 
anism influence particle dynamics only beyond the leading order. The model 
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Figure 1: Diagrammatic representation of the Schwinger-Dyson equation for the scalar theory 
with £i nt = — A(<^T 0) 2 /4. G and Go denote the full and free propagators, respectively, and 
A the four point function. 



whose dynamics we study here is very simple, a complex scalar field with the 
lagrangian 

C = d^d^ - ^m 2 (f>U + h.c. + £ int , (1) 

where for definiteness we shall take the quartic interaction £; nt = — \(<ft <f) 2 / A. 
The mass m = m(x IJ/ ) represents coupling to a classical background field which 
varies in space and time, an important example being the Higgs field conden- 
sate. Such a background implies breakdown of the translational invariance, 
which has as a consequence localization in position space, or equivalently, de- 
localization in momentum space and breakdown of the quasiparticle picture in 
the plasma. This breakdown is quite different from one related to the collec- 
tive self-energy corrections. Indeed, assuming a planar symmetry, m = m{z), 
which models the bubble wall at the electroweak phase transition, we find that 
at the pth order in gradients the momentum splits into p + 2 branches (c/. 
figure 2), while the energy remains conserved. A self-consistent description 
of the problem requires p + 2 distinct distribution functions, or equivalently 
one distribution function and p + 1 functions that measure coherent quantum 
densities on phase space, whose dynamics is described by generalized transport 
equations. In the adiabatic limit (of frequent scattering) the coherent densi- 
ties are suppressed, and one thus recovers the quasiparticle picture with one 
semiclassical transport equation for a distribution function which describes the 
dynamics on a modified shell on phase space. 

2 The Kadanoff-Baym equations 

The Kadanoff-Baym equations (KBE) describe the out-of-equilibrium dynam- 
ics of quantum fields, and can be derived from the Schwinger-Dvson equations 
(SDE) in the Keldysh closed time contour (CTC) formalism □ In figure 1 
we show the SDE for the scalar theory ([!]). In the weakly coupled limit the 
four point function can be approximated by A oc A, so that the SDE for the 
propagator closes. The self-energy E is in this case computed in the Born 
approximation, which is illustrated by the two loop diagram in figure 1. 

Since we are primarily interested in the dynamics of the propagator in 
the presence of a slowly varying background, we write the KBEs in the mixed 
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(Wigner) representation as follows 



cosOjft 2 ±^r}{G r ' a } = 1 (PE) (2) 

-sinO({n 2 }{iG<'>} -{<£<■>}{£*}) 

= icosO ({E>}{G<} - {S<}{G>}) (QBE), (3) 

where the operator 0{f}{g} = \ [dxf ■ dug — dkf ■ dxg] is a generalized Pois- 
son bracket. The PE is the dynamical equation for the retarded and ad- 
vanced propagators, G r ' a , while the QBE (^) is the transport equation for the 
quantum Wigner functions, which in the coordinate representation read 

G>(x,y) = -i( ( f>(x) ( fJ(y)), G<(x, y) = -i(<fl . (4) 

Various quantities in equations (0) and (||) are defined as follows 

Q 2 = Gq 1 -Z R = k 2 -m 2 - E R 
G r ' a = G R t iA, YT' a = Tir t iu;T, (5) 

where and T denote the real part of the self-energy and the damping rate, 
respectively, A denotes the spectral function defined by 

^(C-G a )^(G>-G<), (6) 

while Gr is the real part of the propagator, which is related to the spectral 
function by the following spectral integral Gr = (1/tt)V J dcjfA(oj') / (to — ui'). 
A similar spectral relation holds for Hr and ojT = j(E > — E < )/2. Finally, the 
weak coupling expansion ensures that the self-energies Y* <,> may be computed 
in the Born approximation. 

The mixed representation is defined by the Wigner transform, which is 
the Fourier transform with respect to the relative coordinate r = x — y. For 
example, for the propagator we have 

G{k;X)= J d 4 re lk r G(X + r/2,X - r/2), (7) 

where X is the average coordinate defined by X — (x + y)/2. 

Equations (|^) and (^) represent an accurate description of the plasma 
dynamics background is slowly varying, i. e. when the gradient approximation 
applies. Formally, this is the case when ||0|| <C 1. Physically, this condition 
holds when the de Broglie wave length of excitations is small in comparison to 
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the characteristic length scale of the background variation. For the electroweak 
phase boundary the characteristic scale is the bubble wall thickness L ~ (10 — 
20)/TB, whereas for plasma excitations dk ~ so that ||0|| ~ ||<3x0fe|| ~ 
1/kL, implying that the gradient approximation is accurate when k ^> 
which is close to thermal equilibrium satisfied by most of plasma excitations. 



3 Propagator equation 

We now consider the propagator equation (^|). For simplicity we assume a 
planar symmetry such that to — m(x), and set £_r = 0. For the electroweak 
phase transition this case corresponds to a planar bubble wfill in the wall frame. 
To second order in gradients the propagator is solved bycl 

r r,a^ r 1 , 1 K)" 1 2fc 2 K)" + [(m 2 )f 

P- 2 k 2 — k 2 9 (k 2 - O 0.2 _ JU2\4 ' 

where k% = oj 2 — fcj? — m 2 (x) denotes the unshifted pole, (to 2 )" = d 2 m 2 , and 
(to 2 )' = d x m 2 . Note that at second order in gradients there is no pole shift, but 
instead multiple poles emerge. Now let us make the following decomposition 
of the quantum Wigner functions 

G < = -2iAn, G > = -2iA(n + l), (9) 

where n is the generalized distribution function. Recall that the spectral func- 
tion A (||) defines the spectrum (poles) of the propagator. For example, at 
leading order in gradients the spectral function becomes a sharp on-shell pro- 
jector: A — > Ao = (7r/2fc ) [5(k x — ko) + 5(k x + ko)]. To understand how is 
this projective behaviour modified by the presence of the derivative correc- 
tions in Eq. (0), we study the following spectral integral 

2 f°° P+1 

i p [t] = - / dk x A p r -> ^ Cl r«(fc ). (io) 

Note that the series terminates at the p + 1 st order in derivatives, where p is 
the highest order gradient in the propagator G p . For example, when p = 2, 

XM = ^ + c 2 (dln)(k ) + c 3 (d 3 k n)(k ), (11) 

where k x = k sc defines the semiclassical dispersion relation 

1 (to 2 )" 5 [(to 2 )'1 2 , N 

fc sc = fc + - v ' +— - — (12 
sc u ^ g kl 32 kl y ' 
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This agrees with the dispersion relation obtained by the standard WKB method 
□. The additional higher order derivative terms in Eq. (^IJ) signify the break- 
down of the quasiparticle picture and emergence of quantum correlations, and 
they are direct consequence of the multipole structure of G p (||) . 

4 Transport equations 

We now discuss the implications of the quasiparticle picture breakdown on 
transport. To this end we simplify the QBE (||) bw keeping only the higher 
order derivative term contributing to the flow termH to obtain 

-(0-0 3 /3){^}{ J G<} = \ ({£>}{G<} - {£<}{G>})+o(0 2 ,r). (13) 

In order to make use of the decomposition (^) and the properties of the spectral 
function A, we shall perform spectral integrals weighted by powers of the 
momentum over Eq. ([l3"|). Upon defining the moment functions 

2 f°° 

fl = - dk x k l x A p = 2 n, (1 = 0,1,2,..), (14) 

* J-oo 

these integrals result in the following set of coupled transport equations 

u>d t f t + hd x f l+1 - l -{d x kl)h. x - l SLiM^A {d i k i )fl _ 3 = coil;, (15) 

where Coll; is the collision term, which approximately has the following form 

Coii ; = -r>/, + r<(/, + kl- l m). (16) 

Here T <,> denote the damping rates computed in the Born approximation i 
and, for example, ko = fc sc . To any finite order p in the gradient expansion, 
there are only p + 2 independent moments and in particular at second order in 
gradients one can derive the constraint: f& — 4/co/3 + 6fco/2 — 4feg/i + ^o/o = 0? 
so that Eqs. ( |l5| ) close for the first four moments, I = 0, 1, 2, 3. 

The following picture has thus emerged. Due to the higher order poles 
in the propagator (JsJ) , the quasiparticle dispersion relation has split into four 
branches defined by 

fc - k - k -!)(*- 2) (m 2 )» Z 2 -5* + 5[K)f 

k x = k t - k a + - ^- + - ^— , (17) 

such that for example ko = fc sc in Eq. ([l2]). This is illustrated in figure 2. Note 
that the splitting diverges in the infrared, where gradient approximation breaks 
down. For a particle with an energy u> incoming onto the bubble wall, the 
momentum splits into four branches, and the full description of the dynamics 
then requires defining one distribution function on each branch. 
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Figure 2: The hypersurfaces k x = fc; on which the moments /; flow according to Eqs. ( |li| - |l6| ), 

5 Quantum coherence and semiclassical limit 

We shall now briefly describe an alternative formulation of the transport equa- 
tions (|l5|-^6|), which is convenient for discussion of the semiclassical limit. Let 
us first define 

/ = fcsc/o, /qci = fi — f 

/qc2 = HI K2 - /qcl ~ /, /qc3 = /3A0 — f qc2 ~ /qcl - /, (18) 

which then implies the following transport equation for /: 

uStf + hAU + / q ci) = -r>/ - r<(/ + 1). (19) 

The distribution function / measures population density of particles on the hy- 
persurface k x = k$ = fc sc , while the densities / qc i measure the coherent quan- 
tum correlations between neighboring hypersurfaces k x = and k x = fc, 
defined in Eq. ( |l"7| ) . Eq. (|l^) is already in the form of the standard Boltzmann 
equation, apart from the presence of an unknown function / qc i. Equations for 
fo and /1 then imply 

u(dt + r qc )/ qc i + dxfqcl + d x K2(fqcl + fqc2) + K 2 d x / qc 2 = Si, (20) 

where the source 

si = (fc sc - K 2 )d x f + (d x kl/2k sc ) - d x n 2 f (21) 

represents coherent mixing of / and / qc i- One can obtain similar equations 
for the coherent quantum densities f qc2 and / qC 3. The coherent densities / qc i 
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are all damped at the rate r qc = (r > — T^/lo, which is an out-of equilibrium 
generalization of the on-shell damping rateQ. Finally, in the frequent scattering 
limit defined as 

T qc L > f (22) 

one can neglect / cq i in Eq. ( |l9| ) and one obtains a semiclassical equation in 
which quasiparticles flow along modified semiclassical trajectoricsd. 



6 Conclusions and outlook 

We have presented a self-consistent derivation of transport equations for scalar 
theory beyond the leading order in gradients. We found that as a consequence 
of localization in space, the quasiparticle picture of transport breaks down. We 
have shown how to self-consistently reformulate transport theory by including 
information about the varying background to a finite order in gradients. Our 
result are the transport equations ( |l5[jl^ ), which at second order in gradi- 
ents describe the dynamics on four momentum branches in phase space. We 
have also briefly mentioned that one recovers the semiclassical picture in the 
(adiabatic) limit of frequent scattering. 

The question which we are addressing now is how to generalize our results 
to mixing scalar fields and massive chiral fermions, which are the cases relevant 
for baryogenesis studies. 
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